An analytical model of Faraday rotation in hot alkali metal vapours 
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We report a thorough investigation into the absorptive and dispersive properties of hot caesium 
vapour, culminating in the development of a simple analytical model for off-resonant Faraday rota- 
tion. The model, applicable to all hot alkali metal vapours, is seen to predict the rotation observed 
in caesium, at temperatures as high as 115 °C, to within 1 % accuracy for probe light detuned by 
greater than 2 GHz from the D2 lines. We also demonstrate the existence of a weak probe intensity 
limit, below which the effect of hyperfine pumping is negligible. Following the identification of this 
regime we validate a more comprehensive model for the absorption and dispersion in the vicinity 
of the D2 lines, implemented in the form of a computer code. We demonstrate the ability of this 
model to predict Doppler-broadened spectra to within 0.5 % rms deviation for temperatures up to 
50 °C. 



I. INTRODUCTION 

The ability to control absorption and dispersion in hot 
alkali metal vapours is of great interest in a number of 
areas of atomic and molecular physics. The field of slow- 
light, in particular, is very dependent on this ability 1- 
y]. The control of these properties have other uses in, 
for example, light storage |J|, quantum memory |5|], en- 
tanglement of macroscopic systems Q, and frequency 
up-conversion [7|. In this work we will expand upon 
the flourishing interest in the Faraday effect. This phe- 
nomenon has already been put to many uses including: 
all-optical switching [8|, |9[ , non-invasive atomic probing 
|10| . far off- resonance laser locking JJJ, and realising a 
dichroic beam splitter |12J- A model for absorption and 
dispersion in alkali vapours will provide a tool for the 
prediction and analysis of many of these applications. It 
may also be used to predict the properties of EIT spec- 
tra [13| and slow-light Fourier transform interferometers 
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The Faraday effect is a dispersive phenomenon man- 
ifesting itself in a frequency dependent rotation of the 
plane of linearly polarised light when an axial magnetic 
field is applied to an atomic vapour. This is due to the 
Zeeman shift of spectral lines driven by circularly po- 
larised light. A significant advantage of using disper- 
sive properties in atomic physics is that they dominate 
over absorption in off-resonant schemes [l5|. A model 
has been developed by Siddons et al. for the absorption 
and dispersion in rubidium, which has since successfully 
predicted Doppler-broadened absorption [I61] and disper- 
sion using the Faraday effect [l^- In addition to this 
Marchant et al. have used the Faraday effect in order 
to lock lasers far from atomic resonance [ll[. However 
none of these works characterise the significant tempera- 
ture dependence of Faraday rotation; a property desired 
in many of its applications. 

The motivation of this study is to characterise the 
Faraday effect in caesium in off-resonant schemes. We 
adapt the model produced by Siddons et al. in order 
to predict absorption and dispersion in caesium. This 
model is thoroughly tested through comparison to ex- 



perimental Doppler-broadened absorption spectra. A re- 
sult of this is the identification of a weak-probe regime 
in caesium, in which the effects of hyperfine pumping 
are negligible. Following the successful application of the 
model to the study of absorption, we are confident in the 
model's ability to predict dispersive effects. The model 
for the susceptibility involves Voigt profiles which require 
numerical evaluation using a computer code. However, 
the writing, testing an implementation of such a code is 
complicated and time-consuming. The aim of this work 
is to devise an analytic model for Faraday rotation in off- 
resonant schemes, greatly improving our understanding 
of the phenomenon and aiding applications such as laser 
locking far-detuned from resonance. This paper is struc- 
tured as follows. In section two we present the theory 
behind the existing model. Section three provides details 
on the experimental methods used in the testing process. 
In section four we display our results, discuss their impli- 
cations, and introduce a simple model for off-resonance 
Faraday rotation. Finally, in section five, we summarise 
and conclude. 



II. ABSORPTION AND DISPERSION IN 
ATOMIC VAPOURS 

The transmission of a weak monochromatic light 
source propagating through an atomic vapour is governed 
by the Beer-Lambert law; 



/ = /o cxp [~a {v, T) L] 



(1) 



where / is the transmitted probe intensity for a given 
incident intensity Iq. The absorption coefficient, a, is 
dependent on both the frequency of the light, v, and 
the vapour temperature, T. L is the length of the cell 
containing the vapour. Here we assume that the light 
source is weak enough to not infiuence a. The fraction 
of light transmitted, t, is then 



— = exp I 
Jo 



iL]. 



(2) 



In this paper we are primarily concerned with atom- 
light interactions at frequencies in the vicinity of the D2 



lines for caesium. These lines and the allowed electric 
dipole interactions among them are shown in figure [T] 
Each of the many transitions in this frequency domain 
contribute to a, hence we may express the absorption 
coefficient as a sum 



a i^^, T) = Y1 "j (*"' "^^ ' 
j 



(3) 



where aj {v, T) is the absorption coefficient of transition 
j at frequency v and temperature T. 

In addition to a, each transition also contributes to 
the refractive index, rj, of the vapour. The a^ and a^ 
transitions in an atomic system are driven by left- and 
right-circularly polarised light respectively. When a mag- 
netic field, 5, is applied along the propagation direction, 
the resonant frequencies of the a^ and a^ transitions 
are shifted in opposite directions by the Zeeman effect. 
Consequently the refractive index profiles r?+ and 77^, 
associated with left- and right-circularly polarised light 
respectively, are also shifted in opposite directions. Since 
the linearly polarised light incident on the vapour can be 
considered as a superposition of the two circular polarisa- 
tions, each of these components will therefore propagate 
at different speeds. Upon emerging from the cell, there 
will be a relative phase between the components of 
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(4) 



where A is the wavelength of the probe beam. In the case 
of negligible absorption, the plane of polarisation rotates 

by 



= jL [jf - i-j ) 



(5) 



with respect to the plane of polarisation of the light in- 
cident on the cell. 

The electric susceptibility, x, encapsulates both the 
absorption coefficient and the refractive index, and is a 
complex quantity; 



X = Xr+iXt, 



(6) 



where Xr and Xi ^'^e the real and imaginary parts of the 
susceptibility respectively. The absorption coefficient is 
characterised by 



a (A) = kx^ (A) , 



(7) 



where k is the wavenumber and A is the detuning. 
Throughout this paper, unless specified otherwise, we 
treat zero detuning as the centre of mass frequency of 
the caesium D2 lines. The refractive index is given by 



77(A)-v/xr(A) + l. 



(8) 



In a Doppler-broadened atomic vapour, x for the transi- 
tion i is given as 



X. (A) 
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sj (A) , 



(9) 



where c^ is the transition strength and the A^ is the 
atomic number density [la]- Sj (A) is the lineshape of the 
resonance; a convolution of the Lorentzian atomic line- 
shape and the Gaussian velocity distribution along the 
propagation direction. The total susceptibility is then 
calculated by summing over all transitions. In the first 
half of this work we extend the successful model of x in 
rubidium developed by Siddons et al. 16] by implement- 
ing the same computational methods for atomic caesium. 



III. EXPERIMENTAL METHODS 

A. Experimental set up 

A schematic of the apparatus used is shown in figure [1] 
An extended cavity diode laser at 852 nm was the light 
source used. The combination of a half- wave plate (A/2), 
a polarisation beam splitting cube (PBS) and neutral 
density filters (ND) were used to create linearly polarised 
light of variable intensity. The probe beam then passed 
through a 7.5 cm heated caesium vapour cell based on 
the design of McCarron et al. [13] ■ An external solenoid 
was used to provide both heating, and an axial mag- 
netic field if required. Upon exiting the cell, the probe 
passed through a second half-wave plate and the two or- 
thogonal linearly polarised components of the beam were 
separated using an analysing PBS cube. These compo- 
nents were then focussed onto the ports of a differencing 
photodiode. 
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FIG. 1: (a) A schematic showing the experimental set up. 
Key: half wave plate (A/2), polarising beam splitter (PBS), 
mirror (M), and differencing photodiode (PD). (b) The hy- 
perfine energy states and transitions that make up the D2 
lines, (c) The theoretical D2 transmission in a 7.5 cm cae- 
sium vapour cell at 27 °C (bottom right). The F—4: line is 
shown on the left, and the F—3 line on the right. The con- 
tributions to the absorption of each hyperfine transition are 
shown and labelled. The dashed black curve shows the total 
transmission through the cell. 



The predictions of the computational code were val- 
idated by comparisons to the experimental Doppler- 
broadened transmission of the weak probe propagating 
through hot caesium vapour. These measurements were 
taken under zero magnetic field with the analysing half- 
wave plate set so that all transmitted light was incident 
on a single port of the differencing photodiode. Spectra 
were recorded for both D2 lines at a range of tempera- 
tures using the heating solenoid. 

Faraday signals were obtained by setting the analysing 
half-wave plate so that its optical axes were at 7r/4 to 
those of the analysing PBS cube. The output of the 
differencing photodiode then exhibited zero signal at de- 
tunings where the plane of polarisation was rotated by 
integer multiples of 7r/2. 
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B. Scan Linearisation 

The ECDL used during this experiment incorporated 
current feed forward circuitry. This synchronised mod- 
ulation of the diode current with the piezo allowed scan 
ranges of up to 12 GHz to be obtained. Although the 
laser frequency scans were driven by a triangle wave- 
form, some linearisation was still required. The frequency 
axis of all scans recorded for this work were scaled by 
the use of a Fabry-Perot etalon and saturated absorp- 
tion / hyperfine pumping spectroscopy [18|, [iQj . The dif- 
ferences between the observed and expected detunings 
of the etalon transmission peaks were seen to follow a 
low-order polynomial relation as a function of detuning. 
Subtracting a fit to the residuals from the measured de- 
tuning linearised the scans. The sub-Doppler references 
generated, along with crossover resonances at detunings 
halfway between each pair of hyperfine lines, were at well- 
known frequencies. These features were used to deter- 
mine the absolute detunings across the scan. 



IV. RESULTS AND DISCUSSION 
A. The Importance of a Weak Probe 

The inset to figure [2] shows experimental absorption 
spectra of both the F='i and F=A lines at a tempera- 
ture of 19 °C. These spectra were taken with a probe 
intensity of 0.07 /gat, where the saturation intensity is 
2.7 mWcm"^ 20"]. It is clear to see that, at this inten- 
sity, the predictions of the code underestimate the levels 
of transmission observed in experiment. The F=3 tran- 
sition in particular shows very poor agreement across the 
detunings spanned by the line. 

This occurs even though the probe intensity is below 
the transition's saturation intensity. An underlying as- 
sumption of the model is the cause of this; that the dis- 
tribution of atoms in the hyperfine ground states is un- 
changed by the beam. In reality, it is possible for an atom 



FIG. 2: Linecentre transmission as a function of probe inten- 
sity for the 7^=3 (red circles) and _F=4 (blue triangles) lines. 
The corresponding curves are fits to equation [TO] The hor- 
izontal dashed black lines show the theoretical transmission 
in the weak probe regime for each line. The insets show the 
transmission across each line at 19 °C with a probe intensity 
of 0.07 /aat- The dashed black curves show the theoretical 
transmission. 



excited from one hyperfine ground state, to decay spon- 
taneously into the other ground state, which is known 
as a dark state since is is detuned by 9.2 GHz from the 
resonance being probed. This effect is called hyperfine 
pumping [2l| and is the cause of the reduced absorption 
exhibited in figure [5] 

The high frequency side of the F—'i transition in 
the figure is clearly less strongly affected by this phe- 
nomenon. This is due to the closed transition to F'=5 
dominating the lineshape, as may be seen in figure [1] 
However the weaker closed transition to F'~2 in the 
F=3 line does not have such a significant influence on 
the lineshape. Consequently hyperfine pumping has a 
strong impact over the entire line. 

It has been seen in rubidium, that if the probe beam is 
sufficiently weak, then the time scale of hyperfine pump- 
ing becomes very long and the atoms may traverse the 
beam without being transferred to the dark state [2l| . 
To investigate the effect of hyperfine pumping in the cae- 
sium D2 lines, the linecentre transmission of each line was 
recorded across a large range of probe intensities. The 
results of this study are presented in figure [2] Included 
in this plot are guides to the eye, with the absorption 
coefficient modelled as 



a = Qfwcak ( 1 +/3-p — 
-'sat 



'1/2 



(10) 



where awoak is the theoretical value of a in the absence 
of the pumping. /3 is a fitting parameter characterising 
the change in the saturation intensity due to hyperfine 
pumping. It is immediately apparent from this plot that. 



at intensities below 0.01 /gat, the experimental transmis- 
sion converges to the theoretical value. Despite this it 
is worth noting that this value of the weak-probe limit 
is valid only for the beam used in this experiment, since 
hyperfine pumping will depend upon both the beam size 
and power. However, this is still an important result, 
since we have demonstrated conclusive evidence of the ex- 
istence of a weak-probe regime for the caesium D2 lines. 



B. Comparison of Experiment to Theory 

A series of Doppler-broadened profiles of both D2 lines 
have been recorded in order to test the model against ex- 
periment. These are shown in figure [3] and were taken 
with a probe intensity of 0.001 /sat, which, from figure[2l 
is discernibly within the weak-probe limit. These spectra 
were recorded at a variety of different temperatures using 
the heating solenoid. The spectral lines clearly show the 
expected increase in absorption and width as the tem- 
perature is raised. In addition to the experimental data, 
the theoretical predictions of the model are also shown. 
These are produced by fitting the code output to the 
normalised experimental spectra using the vapour tem- 
perature as a fitting parameter. 

The sub-plots in the figure [3] show the deviations of 
the spectra taken at a temperature of 25 °C from their 
respective fits to theory. It is clear to see that there 
is excellent agreement between theory and experiment 
across both of the D2 lines. The rms deviation is at the 
0.5 % level; a typical level of agreement for the range of 
temperatures used. 

Following minimal adaptation, the code for the elec- 
tric susceptibility of atomic rubidium used by Siddons et 
al. has been shown to successfully predict the Doppler- 
broadened absorption spectra of the caesium D2 lines. 
We have thoroughly tested the code which evidently 
models the quantum system to a good level of accuracy. 
To improve this further, the magnetic octupole interac- 
tion of the caesium D2 lines could be incorporated into 
the model |22| . However this effect has a very minor con- 
tribution to the absorption lineshape and we feel that the 
current level of accuracy is sufficient. To extend the range 
of temperatures in which the model applies, a mechanism 
for coUisional broadenin g wo uld have to be included such 
as that by Weller et al. [23i|. 



C. The Off- Resonance Faraday Effect 

A composite Faraday signal over the D2 lines is shown 
in figure m I^ and ly are the intensities of the two orthog- 
onal linear polarisations analysed by the PBS in figure [l] 
Iq is the intensity of the incident probe beam. For ref- 
erence the absorption spectrum across the D2 lines at 
19 °C is shown in the sub-plot. These plots reveal sev- 
eral interesting features of the Faraday effect. Firstly, it 
is worth noting that the Faraday signal is very sensitive 
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FIG. 3: The left plot shows transmission as a function of 
linear detuning for the ^=4 line at temperatures of 19 °C 
(top curves), 25 °C (middle) and 47 °C (bottom). The right 
hand plot shows transmission for the F—3 line at 21 °C (top), 
25 °C (middle), and 47 °C (bottom). The sub-plots show the 
experimental residuals from the model multiplied by a factor 
of 100 for the 25 °C data. 



to detuning, even in the regions where close to 100 % 
transmission is observed. This in itself makes the Fara- 
day effect extremely useful since signals are produced far 
off-resonance without significant loss of light. Secondly, 
unlike rubidium, the large ground state hyperfine split- 
ting in caesium leads to significant signal oscillation in 
the region between the two lines. This region should 
be ultra-sensitive to the distribution of atoms within the 
ground states. Consequently this frequency domain is 
a candidate for an all-optical switch S, 9']. Finally, one 
notices that the region of zero signal about the linecen- 
tres of the D2 lines does not necessarily correspond to 
100 % absorption. A zero signal may be obtained with 
this optical set up when either the transmission of the 
light driving the a"^ or a^ transition is zero. In both 
cases purely circularly polarised light emerges from the 
cell and will result in zero differencing signal. 

Faraday signals at a variety of different temperatures 
and axial magnetic fields have been recorded. The tem- 
perature dependence of the Faraday effect manifests itself 
in the relationship between the electric susceptibility and 
the atomic number density. Magnetic field dependence 
is due to the splitting of the hyperfine sub-levels via the 
Zeeman effect. The waterfall plots of figure [S] demon- 
strate the evolution of the signals with either of these 
variables. We can immediately see that the vapour tem- 
perature has by far the greater impact on the signals; 
not only do the number of signal oscillations increase, 
but they also are shifted further from resonance. This 
sensitivity is due to the strong temperature dependence 
of the number density and, to a lesser extent, the width 
of the lineshape. Increasing the magnetic field simply 
shifts the edges of the lineshapes, corresponding to the 
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FIG. 4: An example Faraday signal across the entirety of the 
D2 lines at a temperature of 71 °C and a central magnetic 
field of 174±1 G. For comparison a transmission spectrum, 
across the same region, at a temperature of 19 °C is shown in 
the sub-figure. 



translations of the signals seen in figure [S] 



D. Modelling Far Off- Resonance Faraday Rotation 

The Zeeman effect can be included within the model 
described in section two by diagonalisation of the full 
coupling Hamiltonian |l2l |24| | in order to describe the 
Faraday signals. However, it will be shown in this paper 
that a simple model is sufficient to accurately predict the 
extent of Faraday rotation in a far off-resonance regime. 
When greatly detuned from an atomic resonance, disper- 
sion increasingly dominates over absorption. It has been 
shown that, in this off-resonant region, the real part of 
the atomic lineshape factor may be approximated by 



s,(A)«-l/A, 



(11) 



where A is now the detuning from linecentre [15]. Com- 
bining this with equation ini the far-detuned real suscep- 
tibility becomes 



Xr(A) 






(12) 



Since vapour cell number densities are generally very low 
(Xr is typically much less than one) , we can approximate 
the difference between the cr+ and a^ refractive indices 



as 



T]^ 



V 



Xr 



1- 



Xr 



1 



(13) 



Again if we limit ourselves to a regime far from the atomic 
resonances, we may treat all of the constituent hyperfine 
transitions, in the single Doppler-broadened line, as a 
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FIG. 5: Top: Evolution of Faraday signals with a constant 
magnetic field of 117 G with temperatures of 72 °C (red), 
92 °C (blue), and 116 °C (black). Bottom: Evolution of 
signals with a constant temperature of 61 °C and magnetic 
fields of 89 G (red), 145 G (blue), and 210 G (black). 



single transition of strength C^ . Furthermore we assume 
that the Zeeman effect translates into a shift of the re- 
fractive index profiles, ry"*" and r]~ , by a frequency b in 
opposite directions. In this case equation [T51 becomes 



r;+ - 7]- = -C- 



In the limit |A| >> 5; 



T]' -7] 
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(15) 



Combining this with equation [5] the far ofF-resonant Fara- 
day rotation angle is given by 



6 = C^ 



,TTL(fNb 



One of the many apphcations of the Faraday effect is the 
locking of lasers at large detunings from resonance [ll|. 
Consider a set up with which locking points are generated 
at detunings where a rotation in the plane of polarisation 
of n7r/2 is observed, where n is an integer. By this simple 
model, with a constant magnetic field, the locking points 
are expected to occur at detunings 



A„ = 



AN{T) 



(17) 



where ^ is a fitting parameter characterising all the phys- 
ical constants and approximations used, and n is the or- 
der of the locking point. N (T) is a the atomic number 
density as a function of temperature; a well-known func- 
tion for all of the alkali metal atoms 



25| 



In order to test this simple analytic expression for far 
off-resonance Faraday rotation, we recorded the detun- 
ings of the signal zero crossings from the F=A linecen- 
tre as a function of temperature. The detunings of the 
first, third and seventh order zero crossings are plotted 
as a function of temperature in figure [51 where the or- 
der is determined by the integer n in equation 1171 A fit 
of equation [17] to the far-detuned first order zero cross- 
ings is also included. This curve has been scaled to the 
other data sets by dividing by the square root of the ap- 
propriate value of n. The inset shows a Faraday signal 
taken at 115 °C with a central magnetic field of 117±1 G 
and zero crossings corresponding to the data sets on the 
main figure are shown. We can see from the figure that 
at detunings greater than 2 GHz our model works ex- 
cellently, with the data points deviating from the theory 
by no more than 1 %. However, below 1.6 GHz there is 
very poor agreement with percentage deviation exceeding 
30 % in some cases. This behaviour is due to the approx- 
imations used within this theory that are only valid far 
from resonance. In reality the lineshape is much more 
complicated close to resonance than the simple expres- 
sion used in equation I 111 

Another phenomenon observed in figure [5] is that there 
appears to be a temperature dependent threshold detun- 
ing, below which there are no zero crossings. Our model 
fails to predict this since since it does not include the 
dichroism discussed earlier. The orange (light shading) 
region on the plot shows the region in which there is 
less than 5 % a~^ transmission, and in the purple (dark 
shading) region there is more than 97 % transmission. 
These limits are also indicated by dashed lines on the in- 
set. This defines a window in which which zero crossings 
should occur. We can see from the figure that it is indeed 
the case that no zero crossings occur outside this white 
region. 

This simple model works exceptionally well at large de- 
tunings, which is the region it was designed to describe. 
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FIG. 6: The red-detuning from the F=4 line of the 1=' (black 
circles), 3^'^ (red triangles), and T**" (blue diamonds) order 
zero crossings as a function of temperature for a central mag- 
netic field of 117±1 G. The small black points indicate the 
locations of all discernible zero crossings on a 115 °C signal. 
The blue curve is a fit to our model of the 1"* order crossings, 
the other curves are scalings of this fit. The shaded regions 
show the detunings at which there is less than 5 % (orange) 
and more than 97 % (purple) u"^ transmission. The inset 
shows the signal recorded at 115 °C and the locations of the 
various zero crossings. 



This approximation can be used as a quick calibration 
tool for an inexpensive and simple method of locking 
lasers far from resonance; a necessity in optical dipole 
traps. In addition to the reuslts presented here, which 
has focussed on red-detuned light from the F=A line, we 
have verified that this model is successful at other mag- 
netic field strengths for both lines. We have also used the 
model to successfuly predict locking points in the rubid- 
ium D2 lines by fitting to the data presented by Marchant 
et al. [ni. 



V. CONCLUSION 

In summary we have carefully investigated absorption 
and dispersion in hot caesium vapour. We have estab- 
lished the existence of a weak-probe regime in the cae- 
sium D2 lines, in which the process of hyperfine pumping 
is negligible. In light of the identification of this regime, 
the Doppler-broadened D2 absorption spectrum has been 
predicted for temperatures up to 50 °C with a typical rms 
deviation of 0.5 % from experimental data. 

Whilst a code similar to that used in this paper may 
be used to predict Faraday rotation for applications such 
as laser locking, an analytical expression is desirable for 
purposes such as calibration. We have developed a simple 
model that may predict Faraday rotation at large detun- 
ings from resonance. This model has been demonstrated 



to be accurate to within 1 % of experimental data, at de- 
tunings greater than 2 GHz from resonance, for tempera- 
tures up to 115 °C. The major advantage of this model is 
that only a single data set is required for calibration, and 
the other locking points may be extrapolated from a fit to 
the model. Consequently this model has great potential 
as a calibration technique for a simple and inexpensive 
method of laser locking far from resonance. 
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